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Abstract 

We present results of tunneling studies of p-Hgi_ z Cda;Te-oxide-Al structures 
with 0.165 < x < 0.2 in a magnetic field up to 6 T. The tunneling conductivity 
oscillations resulting from the Landau quantization of the energy spectrum 
in the semiconductor volume are investigated. Under an in-plane magnetic 
field the amplitudes of tunneling conductivity maxima connected with the 
tunneling into a and b spin sublevels are found to differ substantionally from 
one another, and the amplitude ratio varying from structure to structure. 
To understand the cause of this behavior, the tunneling conductivity for this 
magnetic field orientation is calculated taking into consideration the multi- 
band energy spectrum. It is shown that the contributions of the different 
spin-sublevels to the tunneling conductivity are dissimilar and the relationship 

between them depends strongly on the value of surface potential. 
PACS number (s): 73.20. At, 73.40Gk, 73.40Qv 
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I. INTRODUCTION 



The electron tunneling in metal - insulator - semiconductor (MIS) structures is a useful 
tool for studying the energy spectrum of two- and three-dimensional states, the interac- 
tion of carriers with the semiconductor-insulator boundary, the interaction with elementary 
excitations such as optical phonons, magnons, etc. [E The tunneling in MIS structures 
based on wide-gap semiconductors, for example GaAs has been studied in detail. [§HT]J 
In these semiconductors the electron effective ^-factor is small and spin-dependent effects 
in tunneling do not manifest themselves. The opposite situation is observed in HgTe-type 
semiconductors with an inverted energy spectrum, where the electrons are the carriers of the 
fourfold degenerate (at k = 0) band. The difference of contributions by different electron 
spin states to various phenomena for such materials was discussed by a number of authors. 
|T2] , p!3|] It has been found experimentally that with the magnetic field parallel to the current, 
the amplitudes of Shubnikov-de-Haas oscillations connected with the different spin states 
differ largely. |12| An analogous difference in tunneling conductivity in structures based on 
inverted-spectrum semiconductors has been observed with magnetic field perpendicular to 



the normal to the plane of a tunneling structure (B _L n). [13| With this magnetic field 
orientation, the energy of the states depends not only on the quantum number, n, and 
the quasimomentum in the direction of the magnetic field, ks, but also oscillation center - 
bounadry distance, xo, and this xq dependence of the energy is essential for oscillation ampli- 
tudes of the tunneling conductivity in the magnetic field. It is well known that for a twofold 
degenerate parabolic band the energy of an electron increases monotonically with decreasing 
Xq. This E-xq relation is the same for both spin states. Consequently, the contributions to 
the tunneling conductivity, that come from tunneling into different spin states, should be 
close in magnitude in MIS structures based on such materials. In terms of the one-band 
energy spectrum model, magnetotunneling under an in-plane magnetic field was considered 



in Refs. |8|-[lT|. In semiconductors with a complex spectrum, the behavior of the different 
spin states near the boundary may differ significantly and it should lead to a distinction in 
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tunneling into different spin states. Thus the behavior of the Landau levels in the vicinity 
of the boundary in the semiconductors with the energy spectrum described by the Dirac 



Hamiltonian was calculated in Ref. ||14|| . It was shown that in a symmetric band-inverted 
junction the Landau eigenenergies vary nonmonotonically with the orbit center - structure 
interface distance and this behavior various for different spin states. The same problem was 
considered also for semiconductors with the inverted spectrum, described by the Luttinger 
Hamiltonian. [15| The behavior of the Landau levels near the boundary is different for dif- 
ferent spin states in this case, too, and this fact was drawn on to explain the significant 
difference of the tunneling conductivity oscillation amplitudes associated with tunneling 
into different spin states in inverted semiconductor p-Hgi.^Cd^Te-insulator-metal tunnel- 
ing structures at B _L n. The fact that the electrons in inverted semiconductors were the 
carriers of the fourfold degenerate r 8 band was essential to the explanation of the observed 
peculiarities. 

Surprisingly, the similar peculiarities of the tunneling conductivity into different spin 
states at B _L n are observed for some tunneling structures based on narrow-gap semicon- 
ductors, where the electrons are the carriers of the twofold degenerate T 6 band. However, 
the ratio of the tunneling conductivity oscillation amplitudes associated with tunneling into 
different spin states varies widely for various tunnel structures. One possibility is that the 
surface potential at the semiconductor-insulator boundary is different in these tunneling 
structures. Such peculiarities of the tunneling into different spin states are the subject of 
the present paper. 

The article is organized as follows. Section [H] presents and analyzes an experimental data 
for tunneling conductivity oscillations with various magnetic field orientations. Section [1 1 1 
is devoted to the calculation of the behavior of the Landau levels near the surface of a Kane 
spectrum semiconductor. The basic equations and formalism are spelled out. Results of the 
calculation of tunneling conductivity oscillations and a comparison with experimental data 
are presented in Section |TV[ 
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II. EXPERIMENTAL RESULTS 



We have investigated the tunneling conductivity = dj/dV and its derivative da^/dV 
as functions of bias V and magnetic field in p-Hg 1 _ a .Cd. c Te-oxide-Al (0.165 < x < 0.21) 
structures with 10 < E g < 100 meV in magnetic fields up to 6 T at a temperature of 4.2 
K. Tunnel structures were fabricated on single-crystal p-Hg 1 „ a ,Cd ;I .Te with uncompensated 
acceptor concentration Na — Np> = (0.8 — 2) x 10 18 cm -3 . The doping level was determined 
from an analysis of galvanomagnetic phenomena in the temperature range 4.2 - 80 K. The 
procedure for the fabrication of tunneling structures was described in our previous papers. 
||16| , [T7|] The resistance of our structures was 0.01 — 1 kOhm. Several tunnel structures 
prepared on each sample were investigated. 

The energy diagram of the tunneling contact is shown in Fig. [l|. The bias shifts the Fermi 
level of the metal relative to that of the semiconductor (Ep) and the tunneling conductivity 
is proportional to the density of states in the semiconductor at energy E = Ep + eV. 
The magnetic field quantizes the energy spectrum of the semiconductor and the tunneling 
conductivity becomes an oscillatory function of magnetic field and bias. Typical oscillation 
pictures of the tunneling conductivity in a magnetic field at a fixed bias for two different 
structures are shown in Fig.|^. It is seen that (i) for both structures the oscillation maxima are 
split, and at B \\ n the amplitudes of the components are comparable in magnitude; (ii) the 
amplitudes of the maxima and magnetic field positions exhibit a rather complicated behavior 
in tilted magnetic field; (iii) at B _L n the ratio of the amplitudes of the components are 
close in magnitude for structure I and they differ significantly for structure II, so that high- 
field components are hard to resolve at this magnetic field orientation. Similar peculiarities 
of oscillation maxima amplitudes are observed in the bias dependences of dj/dV at a fixed 
magnetic field (Fig. |3|). 

A Fourier analysis of the oscillations shows that the cr^ vs B^ 1 curves have two funda- 
mental fields, B\ and B2 = 2B\. The values of fundamental fields do not vary with angle 
between B and n, only the amplitude of Fourier component B 2 decreases when 9 — > 90° so 
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that it practically disappears for structure II. Thus, this behavior of the oscillation provides 
evidence that the tunneling conductivity oscillations for both structures are due to tunneling 
into bulk Landau levels, splitting of the oscillation maxima corresponds to the spin splitting 
of the Landau levels, and the probability of tunneling into spin-up and spin-down states for 
structure II is remarkably different at B _L n. 

The fundamental field B\ at bias V is determined by the value of quasimomentum at 
the energy E = Ep + eV , B{ 1 = 2e/chk 2 , so that measurements at various biases make it 
possible to obtain the dispersion law E(k) for a wide energy range. [l8|jn| For the structures 
investigated, the experimental data are in good agreement with the E{k) curve calculated 
in terms of the Kane model with parameters P = 8.2 x 10~ 8 eV cm, E g = 80 meV, A = oo. 

It is well known that the phenomena occurring in a tilted magnetic field are a very 
complicated problem. We will therefore restrict our attention to the results at B \\ n and 
Bin, namely, the ratio between the oscillation amplitudes for tunneling into different spin 
states and the relationship of this ratio to the structure parameters. 

Let us consider the tunneling current and tunneling conductivity at B \\ n || Ox and 
B J_ n || Ox. For a metal-insulator-semiconductor structure, the Fermi energy of the metal 
is much greater than that of the semiconductor; therefore, upon summation over all metal 
states, the tunneling current for a structure with low barrier transparency can be written in 
the form 

j(V)^Dj2ml=o (!) 

s 

Here D is the barrier transparency, v, \1/ are the electron velocity operator and the wave 
function in the semiconductor, and x = is the coordinate of the boundary. The summation 
runs over all states in the semiconductor, from energy Ep to Ep + eV. 

In a magnetic field, the semiconductor's states are defined by three quantum numbers: 
Landau level number n, oscillator center position, and the quasimomentum component in 
the direction of the magnetic field. At B \\ n || Ox (B = B x ) the quantity (v x ^ n ) 2 depends 
on k x only, and going to integration over energy we obtain 
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j(V)<xDY, ! (vV n ) 2 x=0 dk x 

n •* 

r E F +eV n flh 

= *E/ E , (2) 

/ \ ^0 f \ 2 dkrg . , . 

°d{V) = — CcD {Vty n ) x=0 —\ E =E F+ eV (3) 

In the magnetic field B _L n || Ox. the Landau eigenenergies are a function of the orbit 
center-barrier distance xq and ks, so the square of the velocity at the boundary is also a 
function of x Q and k B , 



j(V)kDJ2 / (v^n)l =Q dx dk B 



/.max 



dE\ 1 



0d(V) oc ^ " («*n)jl=o (^^rj dk B \ E=EF+eV (5) 

Thus, to calculate the tunneling conductivity at B J_ n, one needs to calculate the 
distance dependencies of the velocity and energy near the boundary This problem for a 
semiconductor with a Kane spectrum is considered in the next section. 

III. THE LANDAU LEVELS NEAR THE BOUNDARY AT B_Ln 

To calculate the behavior of the Landau levels near the boundary at B _L n and the tun- 
neling current for this magnetic field orientation, let us consider an insulator-semiconductor 
structure with an abrupt interface (Fig. f|). We assume that the insulator has the same 
energy structure as the semiconductor with the same momentum matrix element P but with 
an energy gap value much greater than that for a semiconductor. The parameters D c , 
and D v are the conduction and valence band offsets, respectively. For such a structure we 
can solve the Schrodinger equation for x < (i.e., for an insulator) and x > (i.e., for a 
semiconductor) independently. Matching these solutions at x = in the required way, we 
find the solution for entire structure. 



We use the Kane energy spectrum model for a semiconductor and an insulator. The 
interaction with remote and spin-orbit split T 7 bands is neglected. We choose the direction 
x to be to normal to the interface and z to be aligned with the magnetic field(5 = (0, 0, B), 
B _L n). The electromagnetic vector potential A is chosen as A = (0,B(x — xq),0). In 
this case the components of the momentum operator are k x = —i-^, k y = —L~ 2 (x — Xq), 



k z = k z = kg. Here L = yhc/eB is the magnetic length. To start with we wish to consider 
this problem for k z = ks = 0. In this case the Schrodinger equation H^> = E^> splits up 
into two independent equations corresponding to the Landau spin sublevels, designated as 
a- and 6- sets. For the a-set these equations read as 



where 



i P ( 8 
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and for the 6-set, 
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H a m a = E^ a , 




H b m b = Em\ 



(6) 
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Here £ = x/L, and £0 = Xo/L are dimensionless coordinates, energies and £7 9 are measured 
from the top of the semiconductor valence band. We seek solutions of (J6|) and (|7|) in the 
form 
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where 



u ± (a, x 2 ) 



-M(a,-,x 2 ) 



„W 1 3 2 , 



and M(a, b, y) is a confluent hypergeometric function. [^T|] Using the properties of the con- 
fluent hypergeometric function, one can write the particular solutions for the a-set 

/ ^ „ \ 



± 



where 



and for the 6-set, 



ELP- 1 ^ (a, (£ - £ ) 2 ) 
ijl u ± ( G _I )( ^_^ o) 2 

iyfiau* (a + ±, (f - £ ) 5 



(9) 
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The solutions for the insulator that correspond to the same energy as the solution for semi- 
conductor are given by (JD, flTOD, where E and E g are changed to E + D v and E g + D c + D v , 
respectively. 

The general solution for x > (R) and x < (L) for the a-set is 
*: L f.r) = ^ a (x) = A"- L #; 
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R,L 



(11) 



s 



The general solution for the 6-set is similar in form. The quantity *f> R has to converge 
at x — > +00 and *f? L has to converge at x — > —00, so we must set A R = B L = 0. Matching 
conditions are obtained by integrating the Schrodinger equation across the boundary. For 
the a-set these conditions are: 



'#(0)1 = tew 



R 



(12a) 



^(0) + ^(0)1 r = U 2 (0) + V3^(0) 



and for the 6-set, 



#(o)i = tew 



(12b) 



(13a) 



V30?(O) +#(0)1 r = [V3#(0) +^ 3 (0) 



R 



(13b) 



These matching conditions give the secular equation for calculating the Xq dependence of 
energy. This problem was solved numerically with parameters corresponding to the samples 
investigated (E g = 80 meV, P = 8.2 x 10~ 8 eV cm) and D c = 2 eV, D v = 1 eV. (Tf 



Such curves for the a- and 6-sets are presented in Fig. [| in dimensionless coordinates 
(x/L,E/hu n ), where uo n = eB/m n c, and m n is the electron effective mass at the conduction 
band bottom. A pronounced difference between the spin sublevels is seen. This is the result 
of the spin-orbit interaction with the boundary. 



IV. DISCUSSION 

It is evident from Eq.fl5|) that to calculate the tunneling conductivity, we need the xq 
dependence of the velocity squared at the boundary. With the definition of the velocity 
operator v = i/h[Hx] we have calculated (v^) 2 for the levels presented in Fig. |5], and as 
an illustration these results for sublevels a\, and b\ are plotted in Fig. |6|. In calculating 
the tunneling conductivity, we have assumed that the Xq dependencies of energy and (v^>) 2 
have the same form at any ks (Fig. [|, g). In Fig. we show the contributions that the 
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tunneling into different Landau sublevels makes to the tunneling conductivity as calculated 
from Eq.(|5]). Note that tunneling into 6-levels is more effective in consequence of differences 
in the E vs xq and (v^) 2 vs xq curves for a and b sublevels. The model discussed here 
has only two parameters, D c , and D v , and the results have only a weak dependence on 
these parameters. Finally, summing the contributions to the tunneling conductivity by 
the different Landau sublevels and subtracting the monotone component, we obtain the 
oscillatory part of the tunneling conductivity presented in Fig. |[ For comparison, we show 
in Fig. H also the oscillations of the tunneling conductivity of the same structure at B || n, 
which are calculated in the same way as in Ref. ||17|| . It is seen that at B || n the amplitudes of 



the tunneling maxima connected with tunneling into a and b sublevels are close in magnitude. 

Inspection of the experimental and calculated results (Figs. [|, ||) shows that the relation- 
ship among the oscillation amplitudes at B \\ n agree well with experimental data for both 
structures. Whereas at B _L n the ratio of the a-maximum amplitude to the ^-maximum 
amplitude is appreciably greater than the calculated value for structure I, but much less for 
structure II. One conceivable reason for this discrepancy is the presence of an attractive or 
a repulsive potential in the semiconductor near the semiconductor-insulator interface. In 
order to verify this assumption the above calculation was carried out for structures with a 
model potential which we suppose as a rectangular well (or a barrier) of width d and value 
U. In this case the wave functions in three regions, x<0,0<x<d and x > d, have the 
same form (pTT]) . It is essential that for < x < d both terms in Eq.([TT]) should be taken 
into account. The requirement that the boundary conditions (|1^) and (|l^) be fulfilled at 
x = and at x = d gives a secular equation for calculating the x dependence of energy. 
The results of calculations for the surface well (ell = —70 meV, d = 25 A) and the barrier 



(ell = 70 meV, d = 25 A) are presented in Figs. £J and [10]. The parameters of the potential 
are so adjusted that the well cannot hold the 2D states. It is seen that the effect of the 
surface potential is different for the a- and 6-sublevels. As a result, the ratio between am- 
plitudes of tunneling conductivity maxima associated with tunneling into a- and 6-sublevels 
is strongly dependent on the surface potential. For the repulsive potential the amplitudes 
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of a- and b- maxima can be close in magnitude (Fig. |10|c), as is experimentally observed for 
structure I (Fig. [3]a), while for a sufficiently large attractive potential only b- maxima are 
visible (Fig. |IO|a), as is experimentally observed for structure II (Fig. |3|b). Note that for 
B || n, the calculated and experimentally observed amplitudes of the tunneling conductivity 
maxima associated with tunneling into a- and 6-sublevels are close in magnitude, for both 
the attractive and the repulsive potential. 

Thus, the ratio of amplitudes of maxima due to tunneling into a and b sublevels at B _L n 
depends on the surface potential and provides a way for estimating the value of the surface 
potential under conditions where the surface potential is small and cannot localize the 2D 
states. 

As mentioned above (ii), the amplitudes of the maxima and the magnetic field positions 
in a tilted magnetic field exhibit a rather complicated behavior (Fig. |2|). This is clearly seen 
for maxima observed in a high magnetic field. The positions of the maxima vary nonmono- 
tonically with angle and with ~ 50° their amplitudes decrease drastically. Basically, the 
same behavior has been observed for all the structures investigated. 

In addition, we would like to point out one further peculiarity of tunneling conductiv- 
ity in a tilted magnetic field. It is most pronounced in bias dependencies of da^jdV in 
the structures based on narrower-gap Hgi-^Cd^Te (Fig. |ll]). A fine structure of maxima 
associated with tunneling into 6-sublevels arises with magnetic field orientations close to 
B _L n (85° > O > 50°). All 6-maxima have a similar fine structure. The calculations of the 
energy spectrum and tunneling in a tilted magnetic field are essential to understanding of 
peculiarities observed. To the authors' knowledge, such calculations for structures based on 
Kane semiconductors are not available. 

V. CONCLUSION 

Tunneling conductivity oscillations in a magnetic field have been investigated for struc- 
tures based on narrow-gap rlgx^Cd^Te. These oscillations in the structure investigated are 
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shown to arise from tunneling into bulk Landau levels. Rather remarkably, at B ± n the 
relationship of the amplitudes of the maxima associated with tunneling into different spin 
states varies significantly for the structures investigated, so that only tunneling into one-spin 
sublevels is observed in some structures. 

To understand the causer of this distinction, we have calculated the behavior of the 
Landau levels near the surface and the tunneling conductivity oscillation at B _L n. The 
relationship of the amplitudes of the maxima associated with tunneling into different spin 
states is shown to depend dramatically on the surface potential. Thus the investigations 
of the tunneling conductivity at B _L n provide a way to estimate the value of the surface 
potential when this is small and cannot localize the 2D states. 
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FIGURES 

FIG. 1. Energy diagram of metal - insulator - semiconductor structure. E c and E v are the 
energies of the edges of conduction and valence bands, respectively. 

FIG. 2. Oscillations of tunneling conductivity in a magnetic field at V = 130 mV for various 
magnetic field orientations for structures I (a) and II (b). 

FIG. 3. The bias dependences of tunneling conductivity in magnetic field 4 T for structures I 
(a) and II (b) for two orientations of magnetic field. 

FIG. 4. The model of the energy structure of insulator - semiconductor boundary used for 
calculations. 

FIG. 5. The xo dependences of the energy of the Landau sublevels in vicinity of the boundary. 
The parameters used are given in the text, B = 4T. 

FIG. 6. The xo dependences of (v^>) 2 for a\ and b\ Landau sublevels calculated with the same 
parameters as in Fig. |, B = 4 T. 

FIG. 7. The contributions to the tunneling conductivity of various Landau levels, B = 4 T. 

FIG. 8. The oscillatory part of the tunneling conductivity calculated for two orientations of 
magnetic field, B = 4T. 

FIG. 9. The xo dependence of the energy of a\ and b\ Landau sublevels calculated for three 
values of the surface potential. 

FIG. 10. Oscillatory part of the tunneling conductivity for three values of the surface potential. 

FIG. 11. The fine structure of the tunneling conductivity oscillations in tilted magnetic field 
for tunnel structure based on p-Hgi-^Cd^Te with E g = 20 meV. B = 5.5 T 
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